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FOREWORD 

This document reports the results of a four-month study conducted 
by personnel of the Lockheed-Huntsville Research &. Engineering Center 
under Contract NAS8- 31800 and entitled "Experiment 74-21 Space Pro- 
cessing Rocket Program" for NASA- Marshall Space Flight Center. The 
stated contract requirements are: 


A. Provide analytical support for fluid 
flow phenomena in Sounding Rocket 
Experiment 74-21. 

B. Evaluate the extent of fluid flow in 
Experiment 74-21 under ground- 
based and sounding rocket flight 
conditions. 


The NASA contract monitor for this investigation was Ms. Carolyn 
Griner of the MSFC Materials and Processes Laboratory (M&PL). Dr. 
May Helen Johnson of MSFC/M&PL served as Principal Investigator for 
Experiment 74-21. 
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1. INTRODUCTION AND SUMMARY 

The Dendrite Remelting Rocket Experiment (Experiment 74-21) was 
performed aboard a Black Brandt VC Sounding Rocket during a period 
which gravity levels of approximately 10 g prevailed. The experiment 
consisted of cooling an aqueous ammonium chloride solution in a manner 
such that crystallization of ammonium chloride crystals proceeded through- 
out a three minute period of "zero-g." The crystallization process during 
flight was recorded on 35 mm panatomic-x film. A number of ground crystal- 
lizations were similarly recorded for comparison purposes. 

The purpose of the present study was to assess the convective and ther- 
mal conditions in aqueous and metallic liquid systems under conditions of the 
flight experiment to help establish the relevance of the rocket experiment to 
metals casting phenomena. The results of the study indicate that aqueous or 

metallic convective velocities in the Dendrite Remelting Rocket Experiment 

-4-5 

cell are of insignificant magnitudes at the 10 to 10 g levels of the experi- 
ment. The crystallization phenomena observed in the Rocket Experiment, 
therefore, may be indicative of how metals will solidify in low-g. The influence 
of possibly differing thermal fields, however, remains to be assessed. The 
Rocket Experiment may also be relevant to how metals solidify on the ground 
at temperature differences and in cell configurations such that ihe flow veloc- 
ities are not high enough to break or bend delicate dendrite arms. Again, 
however, the influence of the nature of the thermal fields must be assessed. 

The possible influence of a dissimilar solute boundary layer also must be 
assessed. Computer calculations of metallic convection at one-g level, while 
possible, involve complications that were beyond the scope of the present study 
to address. The complications are discussed in the Convection Analysis Section. 
Considerations relevant to identifying the conditions whereby a given fluid flow 
pattern of a given velocity and temperature distribution in ammonium chloride 
could be matched in a molten metal system are presented in an appendix. 
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2. EXPERIMENT DESCRIPTION 

The experiment description given here is brief and limited only to 
those aspects that are of interest in the present study. The experiment 
cell had the geometry and dimensions shown in Fig. 1. 


1 mm Thick 
Plexiglas 




+y 


Fig. 1 - Experiment Cell Geometry, Dimensions and Gravity 
Coordinates During Flight 

The coordinates shown in Fig. 1 are the designated gravity vectors during the 
rocket flight. A record of the gravity levels in the designated directions 

during the flight are shown in Figs. 2, 3 and 4. Note that y direction ex- 

- 5 

perienced the highest levels — about 2.5 x 10 g. 

The experiment cell was cooled by thermoelectric devices on two 
sides and the crystal growth viewed through two other sides. The cooling 
and viewing arrangements are shown in Fig. 5. 


2 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 



-Axis Acceleration 


s s 

~ I 

o u 

4 ) 

<M U3 



Fig. 2 - Gravity in X -Direction During Flight Test 



200 



Seconds 

Gravity in Y-Direction During Flight Test 





- Gravity in Z-Direction During Flight Test 





Fig. 5 - Cooling and Viewing Arrangements of Experiment 

Two thermistors, and , were located at the two cooled walls 

and one at the bottom of the cell as shown in Fig. 5. Thermistor Tg was 
separated from the cell wall by 2 .\ mm of plexiglas. The temperature data 
from thermistor E was not utilized in this study hut is included for com- 
pleteness. Assumptions used in this study were that the cooled walls were 
cooled uniformly and the thermistor data accurately reflects the outside 
wall temperature. 

Ground tests were conducted with the cell in two positions with respect 
to the earth gravity vector. The two ground positions are shown in Fig. 6. 

The flight and ground test temperature data are given in Table 1. 

Figure 7 shows how the external temperature, T^, varied in the flight 
test and in the ground-vertical test as a function of time. It can be seen 
that although the rate of temperature change was about the same in the two 
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Ground Based Test Vertical 


Ground Based Test Horizontal 


Light 


Light 


or 



J' 


Fig. 6 - Orientations of Te3t Cell During Ground Tests 
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Table 1 

OBSERVED TEMPERATURES DURING FLIGHT AND GROUND TESTS 


Flight Case 


Time (sec) 

(t = 0, time cooling began 
or 60 sec after liftoff) 


0 

5 

10 
15 
20 
25 
30 
35 
40 
45 
c 0 
'55 
60 
65 
70 
75 
80 
85 
90 
95 
100 
105 
110 
1 15 
120 
125 
130 
135 
140 
145 
150 
155 
160 
165 
170 
175 


Temperature (C) 


t b 

T 1 

T 

‘2 

24.0 

23.8 

23.5 

24.0 

23.5 

23.8 

24.0 

22.9 

22.7 

24.0 

21.3 

19.7 

23.5 

20.5 

17.2 

23.6 

17.4 

15.0 

24.0 1 

15.3 

12.3 

24.0 

15.0 

11.6 

23.6 

12.2 

9.5 

24.1 

11.0 

7.8 

23.8 

9.7 

6.0 

23.5 

8.5 

4.7 

23.3 

6.8 

3.5 

23.3 

6.1 

2.0 

22.9 

5.3 

1.0 

22.7 

3.7 


22.6 

2.7 

I"-' 

22.2 

2.6 

■ ESI 

22.0 

2.0 

-1.3 

21.8 

1.8 


21.5 

1.5 


21.1 

0.5 


20.7 

-0.5 


20.4 

-0.5 

-4.0 

20.4 

-0.6 

-4.1 

20-0 

-2.0 

-4.2 

19.9 

-2.2 

.2 

19-0 

-2.3 

-5.5 

18.7 

-2.7 

-5.7 

18.3 

-2.8 

-6.0 

17.7 

-2.8 

-6.0 

17.3 

-4.0 

-6.3 

16.8 

-4.2 

-6.8 

16.5 

-4.2 

-7.0 

16.1 

-4.6 

-7.2 

16.1 

-4.6 

-7.2 
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Temperature (C) 

T 

B 

T 1 

T 

2 

18.6 

18.6 

16.1 

18.6 

7.2 

7.2 

17.8 

0.0 

0.0 

16.7 

-2.2 

-2.2 

15.0 

-3.6 

-4.4 

13.2 

-4.4 

-5.6 

11.4 

-5.0 

-6.1 

9.72 

-5.3 

-6.7 

8.06 

-5.6 

-6.7 

6.94 

5.56 

-5.6 

-5.8 

-7.7 

-7.5 


0 










?mperature l C ) 
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Fig . 7 - Comparison of Time- Temperature Curves for Temperature 
Tj in Flight and Ground- Vertical Tests 
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cases for about 60 sec, thereafter a marked deviation occurred. The higher 
ground temperatures after 60 sec were primarily the result of vigorous con- 
vection in the liquid. The point is discussed in more detail in the Discussion 
.Section . 

The times at which nucleation was observed in the ground and flight 
tests is as follows: 



Inner T^ Wall or T'j 

Liner T~ Wall or T' 
L 2 

Ground Tests 


oi sec 

Flight Test 


63 sec 


109 sec 

109 sec 
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3. THERMAL ANALYSIS 

Temperatures within the cell as the result of conduction only, where 
g = 0, were obtained from computer calculations. The cell was modeled 
according to the schematic shown in Fig. 8. 




Fig. 8 - Computer Model for Thermal Conduction Analysis 


The numbers in the circles are nodes, i.e., the points between which the 
computer calculates the temperature. Calculated time-temperature values 
at the designated nodes are given in Table 2. The thermal diffusivity used 
for Plexiglas was 1.09 x 10 cm /sec and 1.51 x 10 cm /sec for the 
ammonium chloride solution. The conduction temperature profiles in the 
liquid are generally as shown in Fig. 9. Figure 10 shows temperatures at 
selected nodes as a function of time from initiation of cooldown. 

These temperatures are required for the convection analysis (Section 4) 
and the crystallization analysis (Section 5). 
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Table 2 

CALCULATED INTERNAL TEMPERATURES OF SOLUTION ASSUMING CONDUCTION 
AS ONLY TRANSPORT MECHANISM (ALL TEMPERATURES IN CENTIGRADE DEGREES 


Time After 
Cool r g Initiation 
>ec) 

Observed External 
Wall Temperatures 


Nodes 



Calculated inner 
Wall Temp. 

T i 

T 2 

3 

4 

5 

6 

7 

8 

mm 

mm 

0 

23.8 



23.7 

23.7 

23.7 

23.7 

23.7 

23.7 

■in 

5 

23.5 


23.7 

23.7 

23.7 

23.7 

23.7 

23.7 

23.7 

WTwM 

10 

22.9 

1 22.7 

23.4 

23.4 

23.7 

23.7 

23.7 

23.7 

23.4 


15 

21.3 

19.7 

22.6 

22.0 

23.7 

23.7 

23.7 

23.7 

22.7 

BrTW 

20 

20.5 

17.2 

21.5 

19.7 

23.6 

23.7 

23.7 

23.6 

21.7 

20.1 

25 

17.9 

15.0 

20.1 

17.5 

23.6 

23.7 

23.7 

23.5 

20.5 

18.1 

30 

15.3 

12.3 

18.0 

15.3 

23.5 

23.7 

23.7 

23.4 

18.5 

16.1 

35 

15.0 

11.6 

16.3 

13.5 

23.4 

23.7 

23.7 

23.2 

17.0 

14.4 

40 

12.2 

9.0 

15.0 

12.0 

23.2 

23.7 

23.7 

23.0 

15.8 

13.1 

45 

11.0 

7.8 

13.2 

10.2 

23.0 

23.7 

23.7 

22.8 

14.1 

11.4 

50 

9.7 

8.8 

11.9 

8.8 

22.8 

23.7 

23.7 

22.5 

12.9 

10.1 

55 

8.5 


10.7 

■SI 

22.6 

23.7 

23.7 

22.2 

11.8 

8.8 

60 

6.8 

IS 

9.4 


22.3 

23.6 

23.6 

21.9 

10.7 

7.6 

65 

6.0 


8.2 

■n 

22.1 

23.6 

23.6 

21.6 

9.5 

6.5 

70 

5.3 

1.0 

■ 9 

3.6 

21.8 

23.6 

23.6 

21.2 

8.7 

5.3 

75 

3.7 

9.0 

S 9 

2.9 

21.5 

23.6 

23.6 

20.9 

7.8 

4.6 

80 

2.7 

0.2 


2.5 

21.2 

23.6 

23.5 

20.5 

6.7 

4.2 

85 

2.6 

-0.3 

■ 9 

1.9 

20.9 

23.5 

23.5 

20.2 

6.0 

3.7 

90 

2.0 

-1.3 

4.1 

1.3 

20.5 

23.5 

23.4 

19.8 

5.7 

3.1 

95 

1.8 

-2.2 

3.6 

0.4 

20.2 

23.5 

23.4 

19.4 

5.2 

2.2 

100 

1.5 

-2.3 

3.4 

-0.2 

19.9 

23.4 

23.3 

19.1 

5.0 

1.7 

105 

0.5 

-2.3 

2.9 

-0.4 

19.6 

23.3 

23.3 

18.7 

4.5 

1.5 

110 

-0.5 

-3.4 

2.0 

-0.8 

19.3 

23.3 

23.2 

18.3 

3.7 

1.1 

115 

-0.5 

-4.0 

1.4 

-0.2 

18.9 

23.3 

23.2 

18.0 

3.1 

0.3 

120 

-0.6 

-4.1 

1.2 

-2.0 

18.6 

23.2 

23.1 

17.6 

2.9 

-0.1 

125 

-2.0 


0.7 

-2.2 

18.3 

23.2 

23.1 

17.3 

2.4 

-0.3 

130 

-2.2 


-0.1 

-2.3 

18.0 

23.1 

23.0 

16.9 

1.6 

-0.5 

135 

-2.3 

II £9 

-0.4 

-2.7 

17.6 

23.0 

22.9 

16.6 

1.3 

-0.9 

140 

-2.7 

Hi m 

-0.6 

-3.5 

17.3 

23.0 

22.8 

16.2 

1.1 

-1.6 

145 

-2.8 

-6.0 

-0.9 

-3.8 

17.0 

22.9 

22.8 

15.9 

0.8 

-1.9 

150 

-2.8 

-6.0 

-1.0 


16.7 

22.8 

22.7 

15.5 

0.7 

-2.1 

155 

-4.0 

-6.3 

-1.5 

Hi 

16.4 

22.8 

22.6 

15.7 

0.3 

-2.3 

160 

-4.2 

-6.8 

-2.1 

-4.5 

16.0 

22.7 

22.5 

14.9 

-0.4 

-2.7 

165 

-4.2 

-7.0 

-2.4 

-4.9 

15.7 

22.6 

22.4 

14.5 

-0.6 

-3.0 

170 

-4.6 

-7.2 

-2.6 

-5.1 

15.4 

22.5 

22.3 

14.2 

-0.8 

-3.3 

175 

-4.6 

-7.2 

-2.8 

■ 9 

15.1 

22.4 

22.2 

13.9 

-1.1 

-3.4 

180 

- :.7 

-7.2 

-2.9 

1 9 

14.8 

22.3 

22.1 

13.6 

-1.2 

-3.5 

185 

-4.7 

-8.0 

-3.0 

9 9 

14.5 

22.2 

22.0 

13.3 

-1.3 

-3.8 

190 

-5.0 

-8.2 

-3.1 

m\ u 

14.2 

22.1 

21.9 

13.0 

-1.4 

KQ 

195 

-5.1 

-8.5 

-3.3 

-6.4 

14.0 

22.0 

21.8 

12.7 

-1.6 

Effl 

200 

-6.6 

-8.5 

-3.7 

-6.6 

13.7 

21.9 

21.7 

12.3 

-2.0 

-4.8 

205 

-5.8 

-8.5 

-4.1 

-6.7 

13.4 

21.8 

21.6 

12.1 

-2.4 

-4.9 

210 

-6.0 

-8.2 

-4.2 

-6.6 

13.1 

21.7 

21.5 

11.8 

-2.5 

-4.8 

215 

-5.9 

-8.5 

-4.3 

-6.6 

12.9 

21.6 

21.4 

11.5 

-2.6 

-5.1 

220 

-5.8 

-8.8 

-4.2 

-6.8 

12.6 

21.5 

21.3 

1 1.2 

-2.6 

-5.3 

225 

-6.3 

-8.8 



12.3 

21.4 

21.1 

11.0 

-2.7 

-5.5 

2 30 

-6.6 

-9.3 

Wtmm 


12.1 

21.3 

21.0 

10.7 

-3.1 

-5.8 

235 

-6.6 

-9.3 

sues 


11.8 

21.1 

20.9 

10.4 

-3.3 

-5.8 

240 

-6.F 

-9.2 

-5.0 


11.6 

21.1 

20.8 

10.2 

-3.4 

-5.9 

24 5 

^ 0 

-9.5 

-5.2 

mVm 

11.3 

20.9 

20.7 

9.9 

-3.6 

-6.1 

250 

- . .0 

-9.5 

-5.4 

-7.8 

11.1 

20.8 

20.5 

9.7 

-3.8 

-6.1 

255 

-7.0 

-9.5 

-5.4 

-7.8 

10.9 

20.1 

20.4 

9.4 

-3.8 

-6.2 

263 

-7.0 

-9.5 

-5.4 

-7.8 

10.6 

20.6 

20.3 

9.2 

-3.9 

-6.2 

265 

-7.0 

-9.5 

-5.5 

-7.9 

10.2 

20.3 

20.0 

8.7 

-4.0 

-6.3 

270 

-7.0 

-9.5 

-5.5 

-7.9 

10.2 

20.3 

20.0 

8.7 

-4.0 

-6.3 

275 

-7.0 

-9.5 

-5.5 

-7.9 

10.0 

20.2 

19.9 

8.5 

-4.0 

-6.3 

2' 

-6.8 

-9.5 

-5.5 

-7.9 


20.0 

19.'’ 

8.3 

-4.0 

-6.4 

2,60 

-7.3 

-9.5 

-5.5 


IBS9 

19.9 

19.6 

8.1 

-4.1 

-6.4 

290 

-6.8 

-9.5 

-5.6 

-8.0 

ii9 

19.8 

19.5 

mem 


-6.4 

295 

-7.5 

-9.5 

-5.6 

-8.0 

9.1 

19.7 

19.3 

S&a 


-6.5 

300 

-7.3 

-9.2 

-5.9 

-7.9 

8.9 

19.6 

19.2 

mm 

H 

-6.4 
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4. CONVECTION ANALYSIS 


Dimensionless Analysis : An idea of the types of convection involved 
and their magnitudes can be gained from a comparison of various dimen- 
sionless numbers. The Rayleigh numbers. Ra^ and Ra g are the thermal 
and solutal Rayleigh numbers, respectively. They are defined by the following 
equations 


g 0 T AT d : 


g P c AC d" 


where d = distance between hot and cold walls, g = gravity acceleration, 

(3 ^ = thermal expansion coefficient, (3^. = %nh»te expansion coefficient, 

AT - temperature difference between the hot and cold wall, v = kinematic 
viscosity, a = thermal diffusivity, and = chemical diffusion coefficient. 
The Rayleigh numbers are weighed ratios of the buoyancy force to the 
viscous force and are used to characterize convective flows. For example, 
when the Rayleigh number exceeds certain critical values, heat flow departs 
significantly from pure conductive transfer. In cases of heating from below, 
a critical Rayleigh number defines the point where fluid flow commences. 

In side heating cases, fluid flow always accompanies any finite temperature 
gradient so there is no critical Rayleigh number for onset of fluid flow. 

For the present study these numbers are given in Table 3. Properties 
used in the calculations of Table 3 are given in Table 4. 

Table 3 

C OMPARISON OF RAYLEIGH NUMBERS AT ONE-G 

Dimensionless 

System Conditions Number 


System 

27.8 wt% NH 4 Ct- 


80 wt% Pb-Sn 


1 g, AT = 11 C, d = 0.5 cm Ra T = 3.4 x 10 
1 g, AC = 2.4 wt %, d = 0.5 cm Ra g = 4.9 x 10^ 


1 g, AT = 11 C, d = 0.5 cm 


R T = 591 


1 g, AC = 2.4 wt %, d = 0.5 cm Rg = 1.5 x 10 
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The relative magnitude of the convective velocities must be deter- 
mined with detailed computer calculations because no general analytical 
theory is known that can handle fluids of such high fluidity as metals. For 
a very rough estimate of the fluid velocities for the two fluid systems under 
consideration, use of an equation developed by Batchelor for air^ can be 
used. 


Batchelor's equation for the steady-state velocity distribution for air 
maintained between a constant temperature differential is as follows: 


u 



a R t (1 -- 5 ) (1 



( 1 ) 


where 


u = velocity 

r = distance in liquid from sidewalls 
d = distance between the two sidewalls 
a = thermal diffusivity 
R^. = thermal Rayleigh number 

Maximum velocities occur at one-quarter distances from the side walls. 
Very roughly then, maximum velocities can be given by 

u : 0.01 a R t 
max t 


Similarly 


u = 0.01 D )7 R 
max 1 4 
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Utilizing the foregoing relationships the following velocities are obtained 
for the cases given at the beginning of this section. 


System 

Driving Force 

Estimated Maximum 
Velocities, cm/sec 

27.8% NH 4 Cf . H 2 0 

3.39 x 10 4 

0.5 


R g 4.9 x 10° 

0.9 

80% Pb-Sn 

R t 591 

0.7 


Rg 1.5 X 10 7 

4.5 


On the basis of these figures, it would be expected that solutally iriven 
convection would be much stronger than thermally driven convection and 
that metal systems should show higher convective velocities than ammonium 
chloride solutions. The preceding figures are based on steady state theory. 
The thermal conditions during the rocket flight were not steady state. Thus, 
the solutal velocities obtained in the computer studies show higher solutal 
than thermal velocities. The dimensionless numbers that are relevant to 
estimating how fast maximum convective velocities are attained are the 
Prandtl, Pr, and Schmidt, Sc, numbers. These are defined as follows; 


Pr 


V 

a 



Table 4 gives values of these numbers for the systems involved. Low values 
of these numbers mean that the driving forces for convection propagate faster 
through the liquid than the momentum transfer. Thus higher values of velo- 
city are realized sooner. Because values are generally quite low, solu- 
tally generated convective velocities under the conditions of this study would 
generally be more slowly attained than thermally generated convective velo- 
cities. The preceding discussion illustrates the danger of utilizing dimension- 
less numbers meant to characterize steady-state connection to obtain information 
on transient convection. 
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Another dimensionless number that is of interest in the present study 
io the stability number, Sy, which is defined as follows 

„ Ra S 1 P C AC 

ay ' Ra T Le " p T A t 

where Le is the Lewis number (Le = — ). 

12 

to) 

Lindberg and Haberstrah' ' define the following criterion for the onset 
of a thermal- solutal convection called thermal bursts: 

Sy < yTJe 

It might be mentioned that a variety of thermal- solutal types of convection 
are not clearly defined and characterized at present. Fingering, cellular 
and oscillatory types of convection are possible as the result of the coupled 
thermal and solutal driving forces. The stability number, therefore, should 
be viewed only as a diagnostic aid but not as conclusive evidence that 
thermal -solutal convection does or does not exist in any particular situation. 
For the ammonium chloride system of the present study the magnitudes of 
the stability and Lewis numbers are as follows: 

Sy = 1.8 (AT = 11C, AC = 2.4 wt %) 

Le = 81.1 , /Le"= 9 

Thus, on the basis of the analysis of Lindberg and Haberstrah, thermal 

bursts could well be expected in the ammonium chloride system of the 

present study either on the ground or in low-g. Jet -like convections 

have <>een previously noted in ammonium chloride systems and have been 

(3) 

forwarded as the cause of a defect in metal castings called freckles. 


21 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D496847 


/ 4 \ 

A study of thermal-solutal convection in A#-Mg alloys' ' predicts a 
bund of liquid at the front of a solidifying interface where convection can 
occur while the liquid further out is queiscent. Ai-Mg alloys are similar 
to the ammonium chloride system in that the density decreases ahead of 
the interface as the result of the removal of solute. The value of the band- 
width, x, is such that 




Vx 

exp(-^-) 

12 



G 


where k Q is the equilibrium solid-liquid segregation coefficient and V the 
velocity of crystallization. Convection occurs when 





2.8 


where m is the slope of the liquidus line. The topic of thermal-solutal 
convection is obviously extremely important to metals castings. 

Computer Models; With the aid of the Lockheed Convection Analysis 
Program (LCAP), it is possible to obtain two-dimensional data on internal 
flow patterns, velocities and temperature, concentration, and density dis- 
tributions during transient heating or cooling. In the present study LCAP 
was applied to a variety of problems in connection with the Dendrite Re- 
melting Experiment. The various cases calculated are summarized in 

Table 5. The boundary conditions for the cases are given in Table 6. 
-4-5 

Cases for both 10 g and 10 g were both calculated for two reasons. 
First, the g level data while undoubtedly precise may be somewhat in 
error with regard to accuracy. Apparently some problems exist in cali- 
brating accelerometer data at such low levels. Secondly, because of the 

.4 

small time steps in the computer program required at 10 g, only a few 

•4 

seconds of real time data were available from the calculations. The 10 g 
data, however, was sufficient to ascertain the magnitude of the convective 

_5 

velocities. The 10 g calculations covered a substantial portion of the 
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Table 5 

SUMMARY OK CASES CALCULATED BY LCAP CONFIGURATION 


NH^ Cl solution 

g = 1, thermal 

.4 

g = 10 , thermal, solutal 

-5 

g = 10 , thermal, solutal 

80 wt% Pb-Sn 
-4 

g = 10 and 10 , thermal, 

solutal 

20 wt% Pb-In solution 

g = 1, thermal 

-4 -5 

g = 10 and 10 , thermal, 

sclutal 


® *- 0 — *' 


NH^CI solution 

g = 1, thermal 

.4 

g = 10 , thermal, solutal 

-5 

g = 10 , thermal, solutal 


NH.Cl solution 

’ .4 -5 

g = 1, 10 , and 10 , thermal 


In all of the cases the temperature of the right wall is T^ and that of the left 
wall T'j. Corresponding solute concentrations are and C'^. 
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Table 6 

BOUNDARY CONDITIONS FOR CASES GIVEN IN TABLE 4 


NH4CI Solutions at 10 ~^g and 10~^g; The temperature T'j and T^ were 
obtained from Table 2 for the thermal cases. For the solutal cases the C^ and 
C '2 terms were obtained from a solubility curve of ammonium chloride as a 
function of temperature. The actual times and T and C values used are as 
follows: 


Time (sec) 

Temp. (C) 

Weight (%) 

T ’l 

t' 

2 

C ’l 

C 2 


23.7 

23.7 

27.8 

27.8 

30 

18.5 

16.1 

26.8 

26.3 

60 

10.7 

7.6 

25.0 

24.5 

90 

5.7 

3.1 

24.0 

23.4 

120 

2.7 

0.1 

23.4 

22.7 

150 

0.7 

-2.1 

22.8 

22.2 

180 

-1.2 

-3.5 

22.4 

21.9 

210 

-2.5 

-4.8 

22.2 

21.6 

240 

-3.4 

-5.8 

21.9 

21.4 

270 

-4.0 

-6.3 

21.8 

21.2 

300 

-4.5 

-6.4 

21.6 

21.1 


80 wt % Pb-Sn Solution: About the same thermal and solutal differences 
were utilized as were applied to low-g NH^Cf cases. The initial temperature, 
however, was 300 C. The conditions used correspond tc a realistic case of 
Pb-Sn solidification. 

20% Pb-In: Cases with this system were run to gain insight into how 
anotheT metallic system would behave and to provide an easy experimental 
test of some of the conclusions presented here. 


(Continued) 
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Table 6 (Concluded) 


Time (sec) 

C 

Weight (%) 

T ’l 

t’ 

l 2 

c 'i 

C 2 

0 

260 

260 

20.0 

20.0 

5 

243 

243 

16.7 

16.7 

10 

226 

226 

13.3 

13.3 

15 

210 

210 

10.0 

10.0 

20 

193 

193 

6.7 

6.7 

25 

176 

176 

3.3 

3.0 

30 

160 

160 

0 

0 

35 

160 

160 



40 

160 

160 




NH4 Cl Solutions at 1 g: The thermal convections at 1 g proved so 
violent, as is subsequently discussed, that solutal calculations were not 
pursued. The thermal conditions for the two cases run in the A and B 
configurations are as follows: 



A Configuration 

B Configuration 

Time (sec) 

T ’l 


T' 

1 

T 1 

l 2 

0 

21.7 

22.5 

18.6 

16.1 

30 

15.1 

16.1 

10.4 

11.0 


The above values were obtained by adding 3.2 and 3,8 degrees, respectively, 
to the observed Tj and T^ values. The 3.2 and 3.8 values were obtained by 
subtracting the conduction T^ and (Table 2) values from the observed 
Tj and T^ flight values, i.e., the 3.2 and 3.8 values represent the corrections 
to the outer wall temperatures to obtain inner wall temperatures. 
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time the experiment was conducted and so provided information on transient 
behavior. In Table 5 the A configuration corresponds to the experimental 
cell standing upright while the B and C configurations correspond to it 
lying on its side. The pictorial designations ZZ/Z/and— ► mean insulated 
wall and direction of heat or solute flow. The term ^ g indicates the di- 
rection of the gravity vector. The designations thermal and solutal mean 
that the calculated convection either arises only from thermally generated 
buoyancy or only from solutally generated buoyancy. Neither LCAP nor 
any other known convection program at present can handle coupled thermal- 
polutal cases for times longer than a few milliseconds. The problem with 
the coupled calculation is that of the extremely small time steps required 
to perform the calculation. In order to follow convection for any length of 
time in coupled cases, unreasonably long calculation periods are required. 

A great deal of relevant information, however, can still be gained from the 
uncoupled cases. 

-4 -5 

The 10 g and 10 g ammonium chloride and Pb-In cases show maxi- 

-4 -6 

mum velocities of the orders of magnitude of 10' to 10 cm/sec. The 
Pb-Sn velocities are even lower. Figure 11 shows a couple of typical curves 
for an ammonium chloride solution case. In all of the low-g cases run, the 
internal temperature profiles did not deviate to any detectable extent from 
pure conduction temperature profiles. Figure 12 shows a typical computer 
picture of an internal temperature profile at a giver instant. 

The thermal 1 g ammonium chloride shows maximum velocities of 
11.6 cm/sec. In reality the ground cases were turbulent. The Pb-In system 
was run at 1 g. The time step required for a metallic system is so small 
that only a very small time interval can be calculated in any reasonable length 
of computer time. The 1 g metallic calculations by computer, therefore, 
were not pursued further. It may be mentioned, however, that a solution to 
this problem is available — a new computation technique (alternating direction 
implicit, ADI) could be introduced into LCAP. Such an alteration of the pro- 
gram would take at least six months. 
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Fig. 12 - Computer Picture of Internal Temperature Profile. NH^Ci Solution 
in a Configuration at 10“ g at 53.9 sec After Cooling Initiation 
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Another regrettable feature of not being able to utilize the computer 

to map convection in metallic systems at 1 g is that the metallic systems at 
-4 -5 

10 and 10 g show significant convective oscillatory behavior. Marked 
oscillatory behavior at 1 g could be a significant factor in dendrite remelting. 

A sample problem was conducted to illustrate the behaviors of transient 
convection3. Figure 13 shows plots of maximum velocity versus time for a 
thermal and a solutal case in which the driving forces were set by the follow- 
ing condition: 


0 T AT = p c AC 

i.e., the buoyancy forces are the same in the two cases. The temperature 
or concentration was dropped at time t = 0 at one wall to a constant lower 
value. The physical situation, therefore, corresponds to suddenly applying 
a lower, constant temperature to one wall while the other wall is maintained 
constant at its initial value. Convection will begin and build up until a steady 
state is eventually reached. The purpose of the calculation was to see how 
long it would be before steady state is attained. 

LCAP also provides data such as shown in Figs. 14 and 15. 
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Fig. 13 - Sample Calculation Comparing Convection Development 
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5. DISCUSSION 

Utilizing the crystallization of ammonium chloride from aqueous 
solutions as a model for the solidification behavior of metal alloys in any 
given circumstance requires a comparison of the following factors as a 
function of time: 

• Convective flow velocities 

• Convective flow patterns 

• Temperature and concentration distributions 

• Crystallization kinetics 

The convective flow velocities are of interest in ingot crystallization 
because they conceivably could bend or break off delicate dendrite arms or 
cause various segregation phenomena. The question of how high a velocity 
can be tolerated before bending or breakage occurs has been only briefly 
addressed in past studies. Tiller and O'Hara^ estimate that velocities of 

2 / 

10-10 cm/sec are required to elastically buckle a metallic dendrite raft 

whose stem is 50 to 100 times longer than it is thick. Somewhat smaller 

values are expected for deformation by yield. Further information on the 

criterion for dendrite arm breakage, however, was not found. It would not 

-4 

appear likely that the velocities of 10” cm/sec and lower calculated in the 
present study would be sufficient to break off metallic dendrite arms. Velo- 
cities of 10 ^ cm/sec would be too low to break off ammonium chloride arms 
even though ammonium cnloride crystals are more brittle than metallic 
crystals. The problem requires further study, but this study is inclined 
to discount velocities of 10 cm/sec and lower as causing any dendrite 
arm breakage. 

_4 

Convective flow velocities of about 10 cm/sec can cause a variety 
of segregation behaviors^' ^ The 10"^ cm/sec flow referred to in this 
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case is the interdendritic flow in the mushy liquid-solid region caused by 

solidification shrinkage or expansion with some contribution from gravity 

acting on density variations. Vigorous stirring of the bulk liquid, however, 

was found to be relatively ineffective in producing microsegregation in the 

presence of a metallic dendritic array: ’ The bulk convection velocities of 
-4 -5 

10 to 10 cm/sec calculated in the present study, therefore, can be dis- 
counted as having any significant effect on the crystallization processes. 

The reduction of gravity level undoubtedly reduced also the interdendritic 
fluid velocities. 

As discussed in the preceding section, the calculated convective flow 

-4 -5 

velocities were in the range of 10 to 10 cm/sec for the ammonium 

-5 

chloride cases and 10 cm/sec or lower for the metal cases. In reality 
the velocities in the actual experiment were probably even lower. The 
computer analysis does not consider coupling. The result of coupling would 
decrease convective velocities because thermal cooling of the wall would 
drive the convective flow downward, as shown in Fig. 16a, and dendrite 
crystallization, which depletes the solute concentration, would drive the 
convective flow u* ward, as shown in Fig. 16b. The result could vary from 
the situation shown in Fig. 16c to that shown in Fig. I6e. The exact flow 
pattern is irrelevant to the present study because the uncoupled velocities 
were so low. At higher velocities the coupling could be a possible source 
of shearing action. 

The computer calculations do not consider the heat of crystallization. 
Liberation of .substantial heats of crystallization would have the effect of 
reinforcing the upward flowing convective currents. At low-g this would 
not be expected to be a large effect. Substantial heat liberation on crystal- 
lization would be expected to reinforce remelting. However, because no 
remelting was observed in the low-g test, this effect can also be discounted. 

We now turn our attention to theories of the columar to equiaxed tran- 
sition in metal that do noi depend primarily on gravity -driven convection. 
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c. Heat and Solute 
Removal (One 
Possibility) 



d. Heat and Solute 
Removal (Another 
Possibility) 


e. Heat and Solute 
Removal (Still 
Another Possibility) 
(Ref. ) 


Fig. 16 - Schematic of Various Convections Driven by Thermal 
and Solutal Forces (9) 


There are currently two of these, and we label them for discussion purposes 
as the Supercooling theory and the Remelting theory! ^ The Supercooling 
theory holds that the transition occurs when constitutional supercooling occurs 
ahead of the solidifying interface — nucleation occurs in the supercooled re- 
gion. In the Remelting theory, the small dendrite arms arc melted off by 
thermal fluctuations arising either from a coupj.ed interaction between dif- 
fusional mass and heat transport or thermal, convective oscillations. With 
the Supercooling theory convection is predicted to increase nucleation be- 
cause convective mixing causes the liquidus temperature near the interface 
to increase. Convection would obviously also be predicted to increase equi- 
axed grains in the Remelting theory. In the Dendrite Remelting Rocket 
Experiment no new grains were observed ahead of the solidifying interface. 
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Tentatively therefore, the following conclusions are indicated: the Super- 
cooling theory is wrong and the Remelting theory is right with regard to 
convective thermal oscillations but not coupled diffusional transports. 

These conclusions, howevei , must be tested further with regard to what 
the thermal fields are under low-g and one-g conditions. 

Two theories concerning solute macrosegregation are also candidates 

(12) 

for low-g testing. One of these theories' ’ maintains that interdendritic 
flow caused by solidification shrinkage is the cause of certain macrosegre- 
gation phenomena. The other theory' ' holds that another type segregation 
phenomenon is caused by a mechanism similar to that proposed in the Super- 
cooling theory. The results of the Dendrite Remelting Rocket Experiment 
would tend to refute the latter theory. Again, however, a definite conclusion 
awaits an assessment of the thermal fields. 
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6. CONCLUSIONS 

Gravity-driven convective velocities were insignificant magnitudes 
during the Dendrite Remelting Rocket Experiment. 

Gravity -driven convection would also be negligible in metallic systems 
under the same conditions. 

The Dendrite Remelting Rocket Experiment, showing no secondary 
nucleation phenomena, indicates that Supercooling theories and Remelting 
theories (i.e., remelting by coupled diffusional transport) may not be appli- 
cable to ingot solidification. Definite conclusions, however, await a detailed 
analysis of the respective thermal fields involved in ammonium chloride and 
metallic systems. 

More work is needed on how to match selected convection character- 
istics for both metallic systems and aqueous ammonium chic’ i/\» systems. 
This will allow the aqueous solutions to be utilized more accurately as 
models or analogs for metallic systems. 
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Appendix A 

As discussed in the main body of the report, reduction of convective 
velocities is possible by decreasing (scaling) container dimensions, tem- 
perature differences or gravity levels. Convective velocities, however, are 
not the only parameters that need to be considered in scaling problems. 
Convective flow patterns and temperature/concentration distributions are 
also affected by a change of boundary conditions. For the case of ammonia 
chloride systems as analogues for metallic systems, it will probably not be 
possible to match every dependent phenomenon by a single change of boundary 
conditions, i.e., if we match the flow velocities, the flow patterns or the tem- 
perature profiles will probably be different. Insight into the reasons why an 
exact matching is probably not possible can be obtained from a consideration 
of a dimensionless analysis of the fluid flow equations. This is presented in 
the following paragraphs. The difficulties outlined in trying to match con- 
vective behaviors in two different type fluid systems with highly dissimilar 
thermal and transport properties highlights the advantages of having a well 
constructed computer program. With the aid of the computer the effects of 
a variety of boundary conditions can be quickly assessed. The matching or 
similitude problem then becomes a simple parametric study. 

• THERMALLY DRIVEN CONVECTION SIMILITUDE ANALYSIS 

For gaining some insight into scaling of ammonia chloride experiments 
to Pb-Sn inlow-g flight, the following discussion is presented to help identify 
che important scaling parameters. 

For simplicity, let us look at natural convection on a flat plate using an 
incompressible fluid at steady state. The governing equations become: 


A- 1 


<r 
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Mass 


Momentum 


Energy 



(1) 


( 2 ) 


(3) 


Now for two systems (A and B) of incompressible fluid Eq. (1) will be exactly 
satisfied by simple scaling of the significant lengths and velocities. Let us 
now turn our attention to the momentum equation (Eq. (2)). For system A we 
can write Eq. (2) as follows: 
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Now let us consider another system B related to A by the following: 
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We can write the equation of motion for system B in terms of quantities con- 
tained in system A 
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( 6 ) 


Now the momentum equation for system A (Eq. (4)) is identical to the momentum 
equation for system B (Eq. (6)) if the coefficients of the square bracketed terms 
are equal, i.e., 


C C 2 

P v 


= °T C p C g C P 


C „ C 

M v 


(7) 


We can now substitute the physical quantities (Eq. (5)) into Eq. (7) to obtain the 
following, substituting V and L for significant velocity and length respectively: 


Pb V b/ L B . P B g B B (T ' T J e _ |i B V b /l| 

p A v ^ l a " Pa8a a,T - T «’a ’ "a V L i 

If we combine the first and last terms we find that: 


P B V B L B 


P A V A L A 


( 8 ) 


(9) 


which says that the Reynolds number of system A and system B must be 
equivalent. 


Now combining the second and third terms of Eq. (8) we find that 


P B g B ^B " T -^ L 
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( 10 ) 
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Since velocity is a dependent quantity we eliminate and Vg using Eq. (9) 
and write 


P B *B P B (T ‘ T oo’ b L B 


p a *a ^A < t - V L 1 


which states that the Grashof numbers of systems A and B must be equivalent. 
When buoyancy is the only driving force the fluid velocity is determined en- 
tirely by quantities within the Grashof number and Reynolds number equality 
is superfluous. Therefore, Grashof equality establishes dynamic similarity. 


Now looking at the energy equation (Eq. (3)). 


u 1T + v 8t] - Jjfll 
a* + 8yJ k [ 3 2 y J 


and using the same logic just presented we will show that for similarity of 
temperature fields the Prandtl numbers (C^ JU./k) must be equal for dynamic- 
ally similar systems A and B. Equating coefficients of the bracketed terms 
Eq. (3) in the manner presented for the momentum equations we find 


C p C C C T C v C k C T 


Using the relationship of similar Reynolds numbers for dynamically similar 
flows (Eq. (9)) to eliminate and C dependence we find, after substituting 

physical quantities that; 

C C U. 

_ »A^ A 


which says that for similar temperature fields in dynamically similar flow 
fields we require equal Prandtl numbers. 
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Now since the solutions to the governing differential equations are 
identical for systems A and B if we have equivalent Grashof and Prandtl 
numbers, the Nusselt numbers will be equal and that the magnitude of the 
Nusselt number is only a function of these two quantities, i.e., 

Nu = *(Gr) i|i(Pr) . (14) 

Since the dynamics of the flow are determined by buoyancy forces which, in 
turn, are determined by density changes due to temperature gradients the 
amount of heat transfer (Nusselt number) is the driving parameter of the 
free convection problem. 


If we return to Eq. (7) in our discussion 


C D C V 


°T C p C g C 0 


c c 
M— y 


(7) 


we find that the first term is the coefficient to the inertia forces, the second 
to the buoyancy forces and the third to the viscous forces. 


If we have a case where the inertia forces are small when compared 
with the viscous forces and the buoyancy forces, we can make the following 
similarity statement by neglecting the inertia term. 


C C C 
P g P T 


C „ C 

M v 


Now by making the substitutions (Eq. (5)) we find that 


(15) 
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(16) 


which is the Rayleigh number (Pr • Gr). 
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Now in this case, the heat transfer, Nu, can be defined as only a 
function of this parameter since it completely defines the flow Nu = (Gr • Pr). 
It must be kept in mind that this analysis is good for small inertial forces 
when compared with viscous and buoyant forces. 

-4 - 5 

For the liquid metal case (Pr ~ *01) in space (10 to 10 g) the viscous 
forces and buoyant forces are small and one cannot make the assumption that 
inertial forces are much much less than either of these two. Therefore, the 
equation 


Nu = <p( Gr) j|»(Pr) (14) 

must be used to define the net heat transfer to the fluid. 

The significance of this discussion is, that in order to simulate the liquid 
metal case in low g, the net heat transfer relation (Eq. (14)) must be taken into 
account to achieve similarity in the driving mechanisms. Therefore, one must 
match both Prandtl and Gashof numbers to give total similarity. An additional 
complication arises for the two-dimensional case since the aspect ratio L/h 
must also be equivalent for both cases to achieve similarity. 

The discussion previously presented has shown that dynamic similarity 
(equivalent velocity fields) can be achieved by equal Grashof numbers for sys- 
tems A and B. However, one must keep in mind that the Grashof number can 
be thought of as a "local" Grashof number, that is that it contains spatial de- 
pendence throughout the field. Therefore, if we have two systems, A and B, 
the "local" Grashof number for a low Prandtl number liquid, i.e., liquid metal, 
and high Prandtl number liquid will be markedly dissimilar since the thermal 
energy transport mechanisms of the two systems are completely different, con- 
duction and convection, respectively. This is true even though one may compute 
equivalent Grashof numbers using cell wall temperatures. We can now see that 
the important parameters for dynamic similarity for the one -dimensional case 
are both the Prandtl and Grashof numbers. 
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For the two-dimensional case the added complication of the additional 
momentum equation and added terms to the vertical momentum and energy 
equation require that similarity exists between the aspect ratios of system A 
and system B. 

To summarize this discussion we have found the following to be required 
for dynamic similarity of two thermally driven convection cells: 

• For dynamic similarity one must have equivalent "local" Grashof 
numbers. 

• For equivalent "local" Grashof numbers we must maintain similar 
temperature fields. 

• To achieve similar temperature fields we must have equal Prandtl 
numbers in addition to the similar velocity field. 

• For the two-dimensional case we must have equivalent aspect 
ratios. 

Although all of the convective behaviors of a liquid metal flight experi- 
ment probably cannot be matched exactly with a non-metal model, one or 
another of the convection features probably can be matched with suitable 
scalings. 
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